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ABSTRACT 

The first star formation in tlic Universe is expected to take place within small pro- 
togalaxies, in which the gas is cooled by molecular hydrogen. However, if massive 
stars form within these protogalaxies, they may suppress further star formation by 
photodissociating the H2. We examine the importance of this effect by estimating the 
timescale on which significant H2 is destroyed. We show that photodissociation is sig- 
nificant in the least massive protogalaxies, but becomes less so as the protogalactic 
mass increases. We also examine the effects of photodissociation on dense clumps of 
gas within the protogalaxy. We find that while collapse will be inhibited in low density 
clumps, denser ones may survive to form stars. 
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1 INTRODUCTION 

Althougli we can now detect star-forming galaxies at very 
high redshifts, the ionization and ubiquitous metal enrich- 
ment of the IGM at these redshifts suggest that we are not 
yet seeing the very earliest star formation. In order to study 
the formation of the first stars, we are thus forced to rely 
on a purely theoretical approach. Work over the past few 
decades has led to substantial progress in our understand- 
ing of this epoch, and while many of the details are unclear, 
a broad outline of events can be constructed. 

The first step is the gravitational collapse of overdense 
regions in the early universe. In cold dark matter mod- 
els, this happens hierarchically, with objects on the small- 
est scales forming first. In order for gas to collapse along 
with the dark matter, however, gravitational forces must be 
strong enough to overcome the thermal pressure of the gas. 
This occurs for overdensities more massive than the Jeans 
mass, M j, which, prior to reionizat ion, is of the order of 
10'' M0 dHaiman, Abel fc Rees 2000|). 



As the gas collapses, it is heated by adiabatic compres- 
sion and by shocks. In the absence of cooling, this leads to an 
increase in pressure which will eventually halt the collapse. 
Moreover, it prevents fragmentation of the gas. In order for 
star formation to occur, continued collapse and fragmenta- 
tion are required, and thus stars will form only if the gas 
can lose heat effectively. In quantitative terms, it is neces- 
sary that t he cooling timescale b e shorter than the free-fall 
timescale (Rees & Ostriker 1977). 



During the collapse of the first protogalaxies , the gas 
temperature typically reaches only 1000 - 2000 K (Tegmark 



et al. 1997). This is significantly lower than the temperature 



required for the gas to be able to cool via Lyman-a emission. 
In the absence of metals, the only effective low temperature 
coolant is molecular hydrogen. This forms mainly via the 



reactions 

H + e" 
H" +H 



H" -^7 
H2 -f-e" 



although a small amount also forms via 



H + H+ 
H++H 



W+7 



(1) 
(2) 

(3) 
(4) 



While the final H2 fraction is generally small, it is nonethe- 
less able to cool primordial gas effectively, as has been estab- 
lished by a number of authors (see Haiman, Thoul & Loeb 
1996 and references therein). 

The subsequent evolution of the gas and the eventual 
formation of stars are still not well understood, although 
substantial progress in this area has recently been made 
(Abel, Bryan & Norman 2000; Bromm, Coppi & Larson 
1999; Nakamura & Umemura 1999). 

Finally, once star formation begins, it will influence the 
surrounding gas via a number of feedback processes. Two 
processes in particular have been studied: energy input from 
supernovae and UV radiation from stars. 

Since the energy of a single supernova is comparable 
to the binding energy of the first protogalaxies, it is highly 
probable that they will be entirely disrupted by supernovae 
(MacLow & Ferrara 1999; Ferrara & Tolstoy 2000). However, 
as this occurs only once the first massive stars reach the end 
of their life, it is ineffective during the first few million years 
of star formation. Moreover, the effects of supernovae are 
quite local, being restricted to a single protogalaxy and its 
immediate surroundings. 

Ultraviolet radiation, on the other hand, can potentially 
act much faster, and over a much greater volume. It exerts a 
feedback by photodissociating H2 , thereby suppressing cool- 
ing. The global effects of ultraviolet radiation from the first 
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stars have rec ently been studied by Haiman et al. (200C ) and 
Ciardi et al. (2000). They show that a soft ultraviolet back- 
ground builds up prior to reionization, which photodissoci- 
ates H2 in the IGM and in newly collapsing protogalaxies. 

Less work has been done on the local effects of UV ra- 
diation, ie. the ability of the first massive stars to suppress 
H2 cooling within their host proto galax y. This has recently 
been studied by Omukai & Nishi ([1999), who find that this 



local feedback is very effective, with the radiation from a 
single massive star sufficient to dissociate H2 throughout a 
protogalaxy. However, their treatment of the effects of the 
self-shielding of H2 is correct only once the photodissociation 
region created by the star has reached equilibrium. Without 
knowing the timescale on which this occurs, it is impossible 
to assess the true effectiveness of the local feedback from 
UV radiation. 

In this paper we overcome this problem by introducing a 
simple method to estimate the photodissociation timescale. 
We first calculate the rate at which photodissociating pho- 
tons are produced by a primordial star, approximated as a 
blackbody. Next, we calculate how many of these photons il- 
luminate the object of interest, how many are absorbed, and 
finally, how many of these absorptions lead to photodissoci- 
ation. Most of these quantities can be estimated to within 
a factor of a few or better; even in the worst case, we ex- 
pect our final estimate to be accurate to within an order of 
magnitude. 

Applying this method to simple protogalactic mod- 
els, we find that the photodissociation timescale depends 
strongly upon the mass of the protogalaxy. In the least mas- 
sive protogalaxies, equilibrium is reached on a timescale that 
is short compared to the stellar lifetime, while in larger pro- 
togalajcies, the two timescales are comparable. We also ex- 
amine the fate of H2 in dense clumps of gas within the pro- 
togalaxy. Such clumps are the likely precursors of star for- 
mation, but can collapse only if enough H2 survives within 
them to provide efficient cooling. We find that this is only 
the case within the densest clumps. 

We thus partially confirm the results of Omukai & Nishi. 
Ultraviolet feedback is very effective at suppressing cooling 
in small protogalaxies, if the amount of gas in dense clumps 
is small. It is less effective in larger systems, or when much of 
the gas is densely clumped. This suggests that UV feedback 
alone need not imply a low star formation efficiency. 

In section 2 we examine the manner in which UV radia- 
tion leads to the photodissociation of H2, and briefiy discuss 
the approach taken by Omukai & Nishi and the importance 
of a proper treatment of self-shielding. In section 3 we intro- 
duce our approximate model, and in section 4 apply it to our 
protogalactic model. In section 5 we examine the effects of a 
number of complications not included in our simple model, 
and we conclude in section 6. 



2 PHOTODISSOCIATION OF H2 

The first generation of protogalaxies which are able to cool 
effectively via H2 emission have virial temperatures of the 
order of 1000 K, and mean densities of the order of 1 cm""^ 
Under these conditions. 



in the ground state (in either the ortho or para form, as ap- 
propriate). In order to destroy ground-state II2 molecules, 
a substantial energy input is required. At higher tempera- 
tures (T > 6000 K), this could be provided by collisions, 
but at 1000 K collisional dissociation is unimportant. In this 
case, radiative reactions will dominate the destruction of H2. 
There are several possibilities. 

(i) Photoionization: 

H2 + 7 -» H+ -h e" (5) 

( [O'Neil fc Reinhardt 19781 ). This has a threshold of 15.42eV, 
well above the energy required for Hi ionization. It will 
only occur within the H 11 regi on surrounding the radiation 
source. As we show in section 5.3, this region may be quite 
small. 

(ii) Excitation to the vibrational continuum of an excited 
electronic state of H2: 



H2 -f 7 ^ 2H 



(6) 



( [Allison fc Dalgarno 19691 ). This has a threshold of 14.16eV 
for ortho-hydrogen and 14.68eV for para-hydrogen. This too 
is restricted to Hll regions. 

(iii) Direct excitation to the vibrational continuum of the 
ground electronic state of H2. This process is strongly for- 
bidden and proceeds at a negligible rate. 

(iv) Two-step photodissociation (the Solomon process): 



H2 + 7 ^ 



2H 



(7) 

( ^techer fc Williams 1967 ) . The first step involves excitation 
to either the -B^Sj or C^n„ excited electronic states, with 
threshold energies of 11.15eV and 12.26eV respectively. This 
is followed by decay to some vibrational level of the ground 
electronic state. A fraction of these decays will be to the 
vibrational continuum of the ground state, resulting in dis- 
sociation. A number of rotational and vibrational states of 
-B^E+ and C^Hu are available from the ground state by ab- 
sorption below 13.6 eV; the associated spectral lines belong 
to the Lyman and Werner band systems. 

The fate of H2 in cold, neutral gas thus depends upon 
the flux of UV in the Lyman and Werner bands. Unfortu- 
nately, this is difficult to calculate accurately. Direct solution 
of the radiative transfer problem, even if decoupled from the 
hydrodynamics, is computationally expensive. It is thus use- 
ful to explore analytic or semi-analytic alternatives which 
can give reasonable estimates of the effects of the radiation. 

One possible technique is that suggested by Omukai & 
Nishi. They consider a very simplified chemical model, in 
which the H2 abundance is set by the balance between for- 
mation via H~ and destruction by two-step photodissocia- 
tion. The formation rate is limited by the r ate at which H ~ 
forms. This reaction has a rate coefficient (de Jong 1972) 



fcform = 1.0 X 10 Cni S~ 



(8) 



(Tegmark et al. 1997). Under these conditions, the amount 
of H2 in an excited energy state will be very small, and it 
is a good approximation to consider all of the H2 as being 



for gas at a temperature T. If destruction of H~ by other 
processes can be neglected, then this will also be the H2 for- 
mation rate. This approximation is justified as long as the 
fractional ionization is smaller than 4 x 10~^T^/^ (J. Black, 
private communication), and the density is greater than 
0.045(-Flw/10~^^), where _Flw is the average fiux density 
in the Lyman- Werner bands, in units of erg s~^ cm~'^ Hz~^ 
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( iMachacek, Bryan fc Abel 2000[ ). The first of these condi- 
tions is certainly satisfied; the second becomes so only at 
a distance of the order of 100 pc or more from the nearest 
massive star. However, as we see below, the size of the region 
affected by photodissociation in Omukai & Nishi's model is 
greater still, so their approximation is justified. 

The rate coefficient for H2 photo dissociation can be 
written, in the optically thin limit, as (Draine & Bertoldi 
199e|)~| 

(9) 



fc2stop = 1.13 X lO^Fuw s-\ 



With this simplified chemical model, H2 will reach 
chemical equilibrium on a timescale tdis ~ fc2i3tcp- ^^^^ 
significantly shorter than the timescales on which the fiux, 
temperature, density or ionization vary then we can solve 
for the equilibrium fractional abundance of H2 



^form 

XH2 — "J XqTI 
"'2stcp 



8.8 X lO'^'^XcFLwTn, 



(10) 



where is the fractional ionization and n is the total num- 
ber density of the gas, in units of cm~^. 

If self-shielding of the Lyman- Werner bands is unimpor- 
tant then one can define the 'region of infiuence' of a massive 
star to be the volume of gas surrounding it which satisfies 
two conditions; first, that the cooling time, given the equi- 
librium H2 abundance, is less than the free-fall timescale of 
the gas, and second, that this equilibrium abundance can be 
reached in less than the lifetime of the central star. Omukai 
& Nishi consider the particular example of gas illuminated 
by a massive star with average luminosity density in the 
Lyman- Werner bands of Llw = 10^* erg s~^ Hz~^. In this 
case, using their values of temperature, density and ioniza- 
tion, and assuming a uniform distribution of gas, the cooling 
condition implies a region of infiuence of radius 



^cooi = 110 kpc. 



(11) 



while the equilibrium condition implies a radius of 

req = 9.4 kpc. (12) 

Both of these values are significantly larger than the typical 
virial radius of a protogalactic halo. 



200 



M 



106 Mo 



1/3 



pc, 



(13) 



(where we have followed Omukai & Nishi in adopting a gas 
number density n = 1 cm~^). Hence, Omukai & Nishi con- 
clude that a single massive star can effectively suppress star 
formation throughout the protogalaxy. 

Self-shielding becomes important for gas at a radius r 
from the star when the intervening H2 column density, A'^h2 , 
exceeds 10^* cm^^. In this case the situation is more compli- 
cated. The effect of self-shielding is to attenuate the Lyman- 
Werner flux. For the averaged flux this can be modelled by 
a shielding factor /sh, g iven with reasonable accuracy by 
([Praine fc Bertoldi 199^) 



/sh 



UOi 



(14) 



However, unless the H2 abundance has already reached equi- 
librium, then the column density of H2 between the star 



and a given point, and hence the Lyman- Werner flux at 
that point, will change over time. Omukai & Nishi assume 
that the H2 abundance has reached equilibrium and use the 
previously derived equilibrium abundance to calculate the 
distance from the star at which = lO'^* cm~^. Beyond 
this radius, Fuw drops sharply, and they consider this the 
effective boundary of the region of influence of the star. For 
the same example as above, they find that 



0.97 kpc. 



(15) 



This is still comparable to the size of a typical protogalaxy. 
However, Omukai & Nishi do not calculate the time taken to 
reach the equilibrium state. This may be significantly longer 
than in the optically thin case, as in the initially self-shielded 
gas Flw and hence fc2stcp have values much lower than those 
in the equilibrium case. It is not clear that this equilibrium 
can be reached within the lifetime of the star. 

We attempt to remedy this difficulty by proposing a 
method for estimating tdis for optically thick gas, based upon 
a simple calculation of the time taken for the star to emit 
sufficient photons to dissociate the bulk of the H2. We out- 
line our method in the next section. 



ESTIMATING THE DISSOCIATION 
TIMESCALE 



For a cloud containing a mass 
photodissociation timescale as 



Mh, of H2, we define the 



(16) 



where Mh2 is the mass of H2 destroyed per second. In terms 
of the number of H2 molecules destroyed per second, 7Vh2 
this is simply 

Mh2 = -mu2-^n2, (17) 

where is the mass of a hydrogen molecule. Now, 7Vh2 
can be written as the product of four numbers 

A/'h2 = A^dis/inc/abs/dis, (18) 

where A^dis is the number of photons capable of photodissoci- 
ation produced per second, /inc the fraction of these that are 
incident on the cloud, /abs the fraction of incident photons 
that are absorbed in the cloud, and /dis the fraction of ab- 
sorptions that lead to dissociation. To calculate t^ia we thus 
need only calculate these four numbers. To do this precisely 
is difficult, but fortunately we can estimate them relatively 
easily. 



3.1 Niia 

Any photon with sufficient energy can, in principle, cause 
the dissociation of an H2 molecule, although the probability 
of absorption for those photons with frequencies far from 
the centre of any Lyman- Werner line will be small. To cal- 
culate Adis, therefore, we simply calculate the total number 
of photons between 11.15 eV (the threshold for excitation 
to B^T,^ from the ground state) and the Lyman limit at 
13.6 eV. Higher energy photons will be blocked by neutral 
hydrogen, and are important only within the Hll region. 
For a single star, we have 
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■3^ 

^ CD 




Stellar mass (Mg) 

Figure 1. The number of photons with energies between 11.15eV 
and 13.6eV emitted per second by a metal-free star, as a function 
of stellar mass. The luminosities and effective temperatures of the 
stars are taken from Cojazzi et al. (2000). 



NdiB = / -Adv, 



hu 



(19) 



where Li, is the luminosity density of the star, hui = 
11.15eV and hv2 = 13.6eV. For primordial stars, the metal- 
licity is effectively zero and we treat the stars as black bod- 
ies; for frequencies below t he Lyman break, th is should be a 
reasonable approximation ( Cojazzi et al. 200C ) . In this case 



(20) 



for a star of radius _R, , where F„ = ttB^ is the flux at fre- 
quency f, Bv is the Planck function 



2hV^ 



1 



(21) 



and the stellar effective temperature is Teff. We can write 
this in terms of the total luminosity L, as 

L, 



CII 



Nd\a is then given by 



A^dis 



CIT 



dv. 



(22) 



(23) 



We obtain values of and T^s as a function of stellar mass 
from Cojazzi et al. (2000), who fit an isochrone to a col- 
lection of ZAMS metal-free stellar models. It is then simple 
to calculate A'dis as a function of stellar mass. We plot our 
results in figure Ol. 

We do not incorporate the effects of stellar evolution 
into our model, as we expect them to be a relatively small 
correction. 



• »^ /abs 

The fraction of incident photons absorbed in the cloud is ob- 
viously frequency dependent; far more photons are absorbed 
at frequencies corresponding to the centres of the Lyman- 
Werner lines than in the wings. However, rather than tackle 
this frequency dependence directly, we instead choose to ap- 
proximate the effects of absorption in the set of Lyman- 
Werner lines by 



Wu 



(24) 



where VKtot is the total dimensionless equivalent width of the 
set of lines and Wmax = In (^-[^) — 0.2 is the dimensionless 
width of the range of wavelengths that contains all of the 
Lyman- Werner lines. 

The dimensionless equivalent width of a single line i is 



(25) 



where Tv is the optical depth at frequency v, and the integral 
is performed over the line profile. To calculate VFi, we assume 
that the line shape is well fit by a Voigt profile. We then need 
only a few parameters to calculate the equivalent width. 
We take the oscillator strength of the transitio n and the 
frequency at line centre from Abgrall & Roueff (1989) and 
the intrinsic broadening of the upper lines from Abgrall et 
al. (1992). We also need to know the Doppler broadening 



parameter, 



2fcT 



(26) 



where «t is the microturbulent velocity. Clearly, we will not 
in general know the appropriate value of h to use; however, 
it is possible to place some limits on it. Our lower limit, 
femin = 1.3 kms~^, corresponds to pure thermal broadening 
at a temperature of 200 K. This is approximately the low- 
est temperature at which II2 cooling is effective; the amount 
of gas cooler than this is expected to be negligible, unless 
HD cooling should prove to be important. We obtain an ap- 
proximate upper limit by assuming thermal broadening at 
6000 K (a higher temperature implies significant coUisional 
dissociation), plus a similarly-sized contribution from mi- 
croturbulence. This gives 6max — 10 kms^^, and is a strong 
upper limit. Given some value of h within this range, it is 
then straightforward to calculate Wi as a function of N-^.^ 
by means of the accur ate nu merical approximations given in 
Rodgers & Williams (|l974[). 



For a set of lines, one might naively expect that 



3.2 /inc 

The fraction of the emitted flux incident on a region of inter- 
est depends upon the size of the region and its position with 
respect to the star. For example, for a star in the centre of a 
spherically symmetric cloud of gas, clearly /inc = 1. On the 
other hand, for a star illuminating a compact cloud from a 
distance, /inc = f2/47r, where f2 is the solid angle subtended 
by the cloud, as seen from the star. 



Wtot = ^m, (27) 

i 

ie. that the total dimensionless equivalent width is simply 
the sum of the widths of the individual lines. Unfortunately, 
this expression behaves unphysically in the high column den- 
sity limit, as it does not account for the effects of line overlap. 
Clearly Wtot < W^max - the cloud cannot absorb more light 
than is actually there. To correct for the effects of overlap we 
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10'° 10" 10" 10'° 10'" 10" 

Hg column density (cni~^) 



Figure 2. The fraction of radiation in the Lyman- Werner bands 
that is absorbed by the gas, for Doppler parameter b = femin, for 
an ortho- to para- ratio of zero (soUd hne) or three (dashed Hne). 



follow the procedure suggested in Draine & Bertoldi ( 1996 ) 
using a modified equivalent width defined asQ 



ax }\ ■ 



It is clear from this that Wtot 
then have 



Wn- 



as Wu 



fahs 



Wu 



w„ 



= l-exp{-Wtot/Wm 



(28) 
oo. We 



(29) 



Finally, to calculate Wtot correctly, we must know which 
lines to include. As previously noted, we expect almost all of 
the H2 to be in its ground state. The ortho- to para-hydrogen 
ratio is not known a priori, but there is g ood reason to be - 
lieve that para- hydrogen will predominate (Abel et al. 1997). 
Assuming that there is no ortho-hydrogen, it is easy to cal- 
culate Wtot, and hence /abs, as a function of Nu2- We plot 
this in figure ^ for b = bmin and in figure |^ for b — bmax- 
For comparison, in both of these figures we also plot /abs 
obtained with the assumption of an ortho- to para- ratio of 
three to one. 

We see that uncertainties in the ortho- to para- ratio 
will affect our calculated value of /abs for all but the largest 
or smallest column densities. However, the error introduced 
by this is relatively small - the difference between our two 
extreme cases is never more than a factor of two. The uncer- 
tainty in 6 is potentially more important, but its effects are 
restricted to the square-root portion of the curve of growth. 



between A^'h, 



10"* cm"^ and Nn^ = 10*'' cm""". Even in 



this region, the difference is only a factor of five and outside 
of this region it quickly becomes unimportant. 



3.4 /dis 

The fraction of excitations that are followed by dissociation 
depends upon the vibrational state of the excited lev el; va, l- 
ues for individual levels are given in Abgrall et al. ( 1992 ). 



* Note that while Draine Sz Bertoldi also included the effects of 
absorption in the hydrogen Lyman series in their calculation of 
the effects of overlap, it is clear from their figure 6 that these 
generally have a very small effect on Wtot ■ 




10" 10" 10'° 10'° 10" 10'° 10'" 10'" 10'=' 10" 
Hg column density (cm~^) 
Figure 3. The same as figure 2, but with b = 6max 



Consequently, /dis will depend upon the shape of the inci- 
dent spectrum, and hence the column density of H2. Fortu- 
nately, this dependence is not large. Draine & Bertoldi show 
that, for a power law spectrum, the dissociation probability 
varies between 0.1 and 0.2 with a mean value of 0.15. More- 
over, changing the shape of the spectrum changes this value 
by only a few percent. 

The situation is complicated, however, by the photons 
emitted by those H2 molecules which do not photodissociate. 
The majority of decays are initially to highly excited vibra- 
tional states, and the consequent photons have energies be- 
low 11.15 eV and play no further part in photodissociation. 
A fraction, however, decay to low-lying vibrational states 
and produce photons energetic enough to cause photodisso- 
ciation; indeed, some will decay directly back to the ground 
state, producing photons coincident with ground state ab- 
sorption lines. 

The net effect is to scatter some fraction of the incident 
radiation, while at the same time redistributing it in fre- 
quency space. To investigate the importance of this effect, 
we first assume that the dissociation probability after each 
absorption remains 0.15; this need not be the case, but is 
a reasonable first approximation. Next, we divide the emit- 
ted photons into three classes. First, we have those with 
energies less than 11.15 eV; these can simply be ignored. 
Second, we have those produced by decay to an excited 
state, but with energies greater than 11.15 eV. These will 
not coincide with ground-state absorption lines, and thus 
will only be re-absorbed in the limit of high column densi- 
ties, A''h2 3> 10^^ cm~^, when line overlap becomes impor- 
tant. We denote the fraction of decays producing photons 
of this type by plw- This will depend upon the shape of 
the incident spectrum, but as in the case of /dis we expect 
this dependence to be small. A simple unweighted average 
over all of the excited states available from the ground state 
gives Plw ~ 0.16, which we take as a reasonable estimate. 
Finally, some fraction of the photons are produced by direct 
decay to the ground state. Following Draine & Bertoldi we 
denote this fraction by prot- This is sensitive to the ortho- 
to para-hydrogen ratio, but assuming the standard equilib- 
rium value of three-to-one and performing an unweighted 
average over accessible states gives us pict = 0.08. Alterna- 
tively, if we have pure para-hydrogen, prot = 0.04. We adopt 
the former value hereafter. 
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Armed with these values, we can now calculate the size 
of the correction we must make to /dis- The largest effect 
is found in the limit of high column densities. For A'h2 3> 
10^^ cm~'^, line overlap implies that all photons within the 
Lyman- Werner bands will be absorbed. In this case, we can 
write /dis as: 



his = 0.15 



1 — PLW — Prct 



(30) 



With our values of plw and prct this implies that /dis — 0.2. 
Thus the effect will never be large; indeed, since the size of 
the correction is comparable to our uncertainty in /dis we 
are justified in ignoring this effect in the calculations that 
follow. 



3.5 Calculating the dissociation timescale 

Putting together our various estimates, we have: 



idis — 



(31) 



where 

A/'h2 ~ -'Vdis/inc/abs/dis (32) 

Here is the fractional abundance of H2 and M is the 
total mass of the cloud. Combining the two equations, and 
writing M in units of Mq, we have 



tdis ~ 2 X lO'* 



A'dis/inc/abs/dis 

For the example of a 25M0 star, A'^dis 



10^* and 



tdis ^ 20 



XH^M /AdisV 



fincfahsfdis V 10^^ / 



(33) 



(34) 



4 APPLICATIONS 

As an example, we apply our method to a very simple pro- 
togalactic model, consisting of two components: a spheri- 
cally symmetric, diffuse component (the halo), in which are 
embedded clumps of a significantly denser component. Al- 
though recent numerical simulations offer some support for 
this picture, doubts inevitably remain concerning the gen- 
eral applicability of the results. Moreover, observations of 
the local molecular clouds suggest that it may be an oversim- 
plification. Indeed, there is evidence tha t molecular clouds 
have a fractal struct ure ( Elmegreen 1997 ), al though this re- 



mains controversial (Blitz & Williams 1997). Nevertheless, 



our model has the advantage of simplicity, and should be a 
reasonable guide. 

4.1 The spherical halo 

We first consider the simple case of a single star at the cen- 
tre of a spherical halo. As a representative example, we con- 



sider a 25M0 star, which gives Adis 



10-* 



although it 



is straightforward to rescale our results for a different stel- 
lar mass. Now, spherical symmetry implies /inc = 1, so it 
remains only to specify the mass, fractional H2 abundance 



and H2 column density of the halo. From these we can then 
calculate /abs and /dis. 

We can place a number of constraints upon these pa- 
rameters. The fractional abundance is the most straightfor- 
ward. A lower limit upon its value comes from the require- 
ment that there be enough H2 to cool the gas effectively. 
The required fractional abundance has been estimated by 
Tegmark et al. (1997) to be approximately 5 x 10"*. We err 
on the side of caution and adopt a lower limit of 10"**. An 
upper li mit ca n be obtained via a method suggested by Nishi 
& Susa (1999). They note that the maximum abundance of 
H2 that can form in a halo is approximately the amount that 
forms within some critical timescale. This critical timescale 
is temperature dependent, and is the shortest of the recom- 
bination, cooling and dissociation timescales. For low ion- 
ization gas, this typically gives a maximum abundance of 
order 10"^[] 

We can also constrain the halo mass by a cooling ar- 
gument, given the redshift at which it forms. This depends 
upon the cosmological model. A frequently adopted value is 
z = 30, appropriate to the case of a rare halo in an SCDM 
cosmology. However, values in the range 10 < z < 100 are 
plausible (depending upon the rarity of the halo at forma- 
tion) if we do not limit ourselves to a particular cosmological 
model. 

Given the redshift, a lower limit upon the halo mass 
Ad comes from the requirement that the virial temperature 
of the halo be high enough to form sufficient H2 for effec- 
tive cooling. This has also been studied by Tegmark et al. 
They find that the critical mass is redshift dependent, in- 
creasing at low redshifts. In the range that we consider, their 
minimum value is Mmin = IO^Mq. Adopting a gas to dark 
matter ratio /b = 0.1, this corresponds to a minimum gas 
mass of 10* Mq. 

An upper limit on the mass comes from the fact that 
halos with virial temperatures greater than 10** K will cool ef- 
fectively via Lyman-a emission. Although it is unclear what 
role H2 cooling plays in such halos, and whether it is nec- 
essary for star formation, it is often assumed that Lyman-a 
cooling suffices. As a detailed examination of this issue lies 
beyond the scope of this paper, we also make this assump- 
tion, and do not consider larger halos. The mass correspond- 
ing to a virial temperature of lO^K is given by 



Mn,ax = 2.2X 10^h-\l + z) 



-3/2 



Mr. 



(35) 



where lOO/i km s^^ Mpc^^ is the Hubble constant. In the 
range of redshifts we consider, this corresponds to a maxi- 
mum mass of approximately 10* Mq and hence a maximum 
gas mass of 10^ Mq . 

Given the mass and H2 abundance, we next specify the 
H2 column density, , by choosing a halo density profile. 
The simplest choice is a uniform sphere. In this case the 
spherical collapse model suggests that the resulting virial- 
ized halo will have an overdensity of S — IStt^ with respect 
to the cosmological background at the time of collapse. This 
is equivalent to a total rmmber density of 



t This does not apply to very dense gas, where H2 forms predom- 
inantly via three-body processes and the H2 fractional abundance 
is typically near unity. 
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rivii 



: 2.1 X 10" 



' '1 + 2) cm 



0.0125 J 



(36) 



for a neutral gas of primordial composition, and hence an 
H2 number density 



(37) 



Since we also know the mass, it is now trivial to solve for 
the virial radius of the halo, rvir, and hence for A'h2- We 
find that, for a star at the centre of the halo, the H2 column 
density along a radial ray is 



4.7 X 10' XHoM 



n ^,2 \ 2/3 
15^ j.rl/3 ( ^'b"- \ 



0.0125 I 



(1 



.(38) 



Unfortunately, the uniform sphere is an unrealistic 
model; both theory and computer simulation suggest that 
real protogalactic clouds will be centrally concentrated. To 
address this, we co nsider a second model, that o f a truncated 
isothermal sphere ( Shapiro, Iliev fc Raga 199S: ). This has a 
core overdensity S = 1.796 x 10^, a core radius ro — 0.038rvir 
and a truncation radius n — l.lOSrvir (corresponding to 
rt = 29.4ro). Over much of its radius, the profile is approxi- 
mately isothermal. Comparing the two column densities, we 
find that 



TVt = 11.3Nu, 



(39) 



where A'r and A''u are the column densities for the truncated 
isothermal sphere and uniform sphere respectively. The uni- 
form sphere likely underestimates the column density by an 
order of magnitude. However, in the truncated isothermal 
sphere model, a significant fraction of this column density 
comes from the small, dense core. Outside of the core, the 
uniform sphere model may give a better estimate. We there- 
fore continue to consider both models, with the expectation 
that the uniform sphere model will provide us with a strong 
lower limit on tdis, while the truncated isothermal model 
should provide a reasonable upper limit. 

There remain two parameters to specify, related to the 
microphysics of the gas. These are the Doppler parameter b 
and the ortho- to para-hydrogen ratio. Both of these must 
be regarded as rather uncertain. We consider two sets of 
extreme values. For the uniform sphere model, we assume 
b — &min and an ortho- to para- ratio of zero (ie. no ortho- 
hydrogen) , while for the truncated isothermal sphere model, 
we assume b = bmax and an ortho- to para- ratio of three. 
Our two models now give strong upper and lower limits on 

Applying our calculations to the concrete example of a 
halo of gas mass M — 10^ M©, with an H2 abundance of 
5 X 10"'* that formed at z — 30, we find that 



8.3 X 10* yr < tdis < 9.5 x lO'"^ yr. 



(40) 



To put these figures into context, the massive stars we are 
considering will have main sequence lifetimes that are typ- 
ically of the order of a few Myr, so in this particular case 
cooling will cease before the end of the star's life. 

We next proceed to vary the parameters individually, 
to examine their relative importance. We take as a fiducial 
model the example above, vary M, XH2 and z over their 
allowed range, and plot the results in figures ^ to |^ respec- 
tively. Several points are immediately apparent. 

Firstly, the most important parameter for determining 




10° 10° 10' 

Gas mass (Mq) 

Figure 4. Upper and lower limits on the photodissociation 
timescale, plotted as a function of the mass of gas in the halo, 
for halos forming at z = 30 with H2 abundance xjj, = 5 x 10"*. 



tdis is the mass of gas in the halo. Its importance is due to the 
weak dependence of A'^h2 upon mass. As the mass of the halo 
increases, the mass of H2 within it increases proportionately, 
but the column density increases at a much slower rate. This 
means that a large increase in mass leads to only a small 
increase in the amount of radiation absorbed, and hence the 
increase in tdis is large. Moreover, our constraint upon the 
halo mass is much poorer than those upon redshift or H2 
abundance. 

Another conclusion that can be drawn is that changes in 
z (which are equivalent to changing the H2 column density 
while keeping the mass of H2 in the halo fixed) have only a 
small effect. This is because the H2 column densities for the 
halos considered generally correspond to equivalent widths 
in the square-root portion of the curve of growth. Thus, large 
changes in A^H2 lead to only small changes in /abs and tdis- 

It is clear from figure ^ that in the smallest halos pho- 
todissociation equilibrium will be reached quickly, within a 
small fraction of the stellar lifetime. We can therefore apply 
the results of Omukai & Nishi with confidence to such ha- 
los. As the mass of gas in the halo increases, however, tdis 
becomes steadily more important, and the assumption of 
equilibrium breaks down. For example, when the gas mass 
is 10^ M0, the lower limit on tdis is approximately 10® years - 
a significant fraction of the stellar lifetime. We cannot safely 
apply Omukai & Nishi's results to such halos, and they may 
undergo significant cooling right up until the point at which 
the massive star becomes a supernova. 

This situation changes somewhat if we increase A'^dis, 
as would occur if a more massive star formed, or if several 
massive stars formed at approximately the same time. In 
this case our estimated values of tdis would be correspond- 
ingly lower, and cooling would be suppressed much sooner. 
Nevertheless, even in this case we cannot conclude that UV 
feedback necessarily implies a low star formation efficiency, 
as much gas may already have cooled and collapsed into 
dense clumps by the time the first star forms. To properly 
ascertain the effects of the ultraviolet radiation on star for- 
mation efficiency, we must also consider the fate of H2 within 
these dense clumps. 
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Figure 5. As figure y, but examining tlie effects of varying tfie 
redsliift of formation. Tlie gas mass is fixed at IO^Mq, and tfie 
H2 abundance is again = 5 X 10~^. 




SxlO * 5x10 
Hg abundance 

Figure 6. As figure ^, but examining ttie dependence upon H2 
abundance. Tlie gas mass is 10^ Mq and tfie fialos form at a red- 
sliift z = 30. 



4.2 Dense clumps 

The second component of our protogalactic model consists 
of dense clumps of gas embedded in the diffuse halo. The 
free-fall timescale of such clumps, ts, will be considerably 
shorter than that of the halo. In the absence of significant 
pressure support, we would expect the clumps to collapse 
quickly and presumably to form stars on a similar timescale. 
This will only be true so long as the clumps remain able to 
cool efficiently. We assume that clumps will collapse and 
form stars only if they retain their H2 until collapse, ie. only 
if tdis > t«- 

To see whether clumps are capable of preserving their 
H2 for long enough to form stars, we need to calculate both 
of these timescales. The free-fall timescale is straightforward; 
for gas of density p it is simply 



iff 



inS -1/2 

10 n ' yr. 



(41) 



Comparing this with a typical massive star lifetime of a few 
Myr, we see that we need not consider clumps with densities 
lower than lO^cm"^; at such densities, clumps will not have 



collapsed by the time the star becomes a supernova. Their 
subsequent evolution will be shaped by their interaction with 
the supernova remnant, a topic which lies beyond the scope 
of this paper. 

Our calculation of tdis for denser clumps is complicated 
by the fact that at high densities our assumption that all H2 
molecules will be in the ground state is no longer valid. At 



densities above some critical value ricr, coUisional excitation 
and de-excitation dominate, and the populations of excited 
states approach their LTE values. The value of the critical 
density depends upon the particular excited state, as well as 



the temperature, but typically 



10'* cm^ for the lower 



rotat ional states (Le Bourlot, Pineau des Forets & Flower 



199!;). Vibrational states have a much higher critical density 
(ricr ~ 10** cm^"*), and in any case will have a negligible 
population at 1000 K. 

In calculating tdis for dense clumps we therefore assume 
that the rotational states are in local thermodynamic equi- 
librium, with a temperature of 1000 K. The main effect of 
this is to increase /abs (and hence decrease tdis), typically 
by a factor of two to three. 

Next, we note that for dense clumps /inc ^ 1; indi- 
vidual clumps receive only a small fraction of the total flux 
from the star. Thus, despite having a larger value of /abs 
and a much smaller mass of H2 than the typical halo, it 
does not necessarily follow that a clump will have a shorter 
dissociation timescale. 

For a clump of angular size Q., as seen from the star, 
we have /inc = il/47r. In the particular case of a spherical 
clump, in the small angle approximation. 



/inc 



(42) 



where R is the radius of the clump, and D is the distance 
from clump to star. For simplicity we do not consider ellip- 
soidal or filamentary clumps; however, it would be straight- 
forward to extend our analysis to them. 

Modelling the density structure of the clumps is prob- 
lematic; we have little theoretical guidance. Consequently, 
we adopt the somewhat unrealistic model of a uniform den- 
sity sphere. This will minimize tdis and maximize ta for a 
given clump. We can then be confident that if our model 
clumps collapse, so would more realistic clumps. 

For a given clump we specify four parameters: its mass 
(in units of M0), density, H2 abundance and distance from 
the star. From these quantities, we calculate the clump ra- 
dius 



i? = 6.2 X 10"'Ar^-'n 



18 » ri/3„-i/3 



(43) 



and the H2 column density along a radial ray passing 
through the centre of the clump. 



Nho = 1.2 X 10*^a;HoM*''='n^/^ cm 



(44) 



From these we calculate /inc, /abs and /dis and thence tdis 
We find that 



tdis = 20.H.M*/^n^/^O^ (^)"Va-bt/dis^ yr, 



(45) 



where D is in parsecs. 

We immediately see that the dependance of tdis upon 
clump mass is weak, unlike the case of the halo. This is 
due to the scaling of fine with clump mass; since massive 
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clumps are larger, they absorb more of the flux, which, to 
a great extent, offsets the fact that they contain more H2. 
This weak dependence is rather fortunate, as we know very 
little about the mass spectrum of clumps to be found within 
protogalaxies. 

The dependence of tdis upon is stronger, but we do 
not expect the H2 abundance in clumps to differ much from 
that in the halo; consequently, the same tight constraints ap- 
ply, and the overall uncertainty in tdis is small. In very dense 
clumps, as noted previously, this does not apply - three- 
body formation of H2 allows much higher abundances to 
form. However, this process only becomes important above 
a number density of 10^ cm~"^; as we show below, in clumps 
of this density, H2 will in any case survive until collapse. 

The most important parameters for determining the 
fate of a clump are thus its density and its distance from 
the star. For a given mass and H2 abundance, we can set 
tdis equal to and solve for the critical distance Dcrit at 
which H2 survives until collapse; this will be a function of 
clump density. We examine three cases, with different values 
of M and Xii2- Rather than varying these quantities sepa- 
rately, we instead vary the combination Xii2M^^^ , as it is in 
this form that the mass and H2 abundance appear in our 
expressions for the H2 column density and for tdis- We ex- 
amine cases corresponding to values of 10"*, 10~^ and 10~^ 
for x-ii2M^^^; we plot the results in figures ^, ^ and ^respec- 
tively. These values cover the range of clump masses and H2 
abundances that we would reasonably expect to encounter. 

The qualitative behaviour is clear. High density clumps 
easily preserve their H2 until collapse, unless they happen to 
be very close to the star. Lower density clumps, on the other 
hand, are prevented from collapsing at significant distances 
from the star - to put these values for Dcrit into context, the 
virial radius of a 10^ M0 halo collapsing at 2 = 30 is about 
100 pc. 

To quantify the effects of ultraviolet feedback upon star 
formation within clumps is rather harder, as in this case we 
are not interested in the absolute value of Dcrit so much 
as its size relative to the size of the star-forming region. 
Nevertheless, a rough guide to the expected effects would 
come from assuming that stars will form only in clumps 
denser than some threshold, say n = 10® cm~'^, for which 
Dcrit is small compared to the size of the protogalaxy. If we 
increase A^dis (by forming more massive stars, for example), 
then this density threshold will increase. On the other hand, 
we expect our model to give us an upper limit on Dcrit and 
hence on the density threshold; realistic clumps may well 
survive much nearer to massive stars than our model clumps 
do. 



5 COMPLICATIONS 

In order to construct a simple approximation, we have omit- 
ted a number of complicating features. Some of these omis- 
sions call for further comment. These are the heating of the 
gas during photodissociation, the effects of circumstellar ma- 
terial, the contribution to H2 dissociation made by Hll re- 
gions and stellar winds and the possible effects of dust. 




1000 10* 10^ 10^ 10' 10' 

Number density (cm~^) 



Figure 7. The distance at which the dissociation and free-fall 
timescales are equal, plotted as a function of clump density, for a 
clump with XHaAfl/^ = 10~*. 




1000 10* 10^ 10^ 10' 10° 
Number density (cm~^) 

Figure 8. As figure R, but for a clump with xn^M^^'^ = 10~^. 




1000 10* 10^ 10^ 10' 10" 
Number density (cm""^) 

Figure 9. As figure ^, but for a clump with xn^M^^^ = 10~^. 

5.1 Gas heating due to photodissociation 

Rather than simply suppressing cooling, H2 photodissocia- 
tion will in actual fact lead to heating of the gas, as some 
of the energy of the Lyman- Werner photon is transferred 
to the kinetic energy of the hydrogen atoms. However, it is 
easy to show that the the effects of this heating upon our 
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protogalactic model will be small. A si ngle photodissocia- 
tion typically produces 0.4 eV of heat (Black & Dalgarno 
1977), ind thus destruction of all of the H2 within a cloud 



will increase the temperature (in the absence of any other 
heating or cooling) by 

6.4x10-^2 
AT = a;H2 K 

~ 5OOOOXH2 K. (46) 

As our upper limit on x^i^ is 10~^, this corresponds to a 
maximum temperature increase of roughly 50 K. This is too 
small to significantly affect our results, and thus we are jus- 
tified in ignoring this effect. 

5.2 Circumstellar material 

Since we expect stars to form within dense clumps of gas, 
it is reasonable to assume that they will still be surrounded 
by large amounts of molecular gas after they have formed. 
If this is so, then this gas may absorb most of the Lyman- 
Werner photons emitted by the star, in which case dissoci- 
ation in the halo or in other dense clumps will be reduced 
until the bulk of the circumstellar H2 is destroyed. 

If this material is distributed in a spherically symmetric 
fashion around the newly formed star, then it is quite easy 
to estimate its effects. There are two possibilities. If the bulk 
of the gas is of moderate density, then we can simply adapt 
our estimate of tdis for a spherically symmetric halo to the 
case of a far less massive clump. In this case, it is clear that 
the dissociation timescale will be short. On the other hand, 
if the circumstellar gas is dense enough to form H2 via three 
body reactions, then xn^ ~ 1, and the amount of H2 gas can 
be large. However, such a high density also implies /abs — 1, 
which goes some way to compensating for this. We find that 
in this ;ase: 



tdv. 



(47) 



This is short compared to the stellar lifetime for any plau- 
sible value of the clump mass. 

However, both of these cases are rather unrealistic. Not 
only do they neglect the effects of ionizing radiation and 
stellar winds, which can be expected to help clear away any 
surrounding molecular gas, but they also assume spherical 
symmetry. In reality, we would expect the gas surrounding 
newly formed primordial stars to be highly flattened. This 
is necessary for efficient outward transportation of angular 
momentum, in the assumed absence of magnetic fields. In 
this case, we would expect the circumstellar gas to survive 
for a much greater time, as it will receive less of the stellar 
fiux. Nevertheless, it will not seriously affect our estimates 
as it will shield only the small fraction of gas which lies in 
its shadow. 



5.3 Hll regions and stellar winds 

Both ionizing radiation and stellar winds offer plausible al- 
ternative mechanisms by which H2 may be destroyed, and 
which may reduce the importance of H2 photodissociation. 
It is thus important to examine their effects. 

Stellar winds from low metal licity stars have recently 
been studied by Kudritzki (2000). They are found to be 



several orders of magnitude less efficient than winds from 
solar metallicity stars, due to the lack of metal lines with 
which to drive them. It thus seems reasonable to neglect 
their effects. 

Ionizing radiation, on the other hand, is rather more 
important. Ionizing fiux from a massive star will lead to the 
formation of an Hll region surrounding the star, whose size 
and subsequent evolution depend upon the density distri- 
bution near the star. A general treatment of the formation 
and evolution of this H 11 region requires solution of the cou- 
pled radiative transfer and hydrodynamics problem, which 
is clearly impractical. However, we can get some guidance 
as to the expected effects by examining some relevant cases 
where the high degree of symmetry allows for analytical so- 
lution. 

The simplest case is that in which the gas surrounding 
the massive star is spherically symmetric and of uniform 
density. In this case, th e growth of the H 11 region can be 
divided into two phases ( Yorke 1986| ). During the first phase, 
an R-type ionization front sweeps into the surrounding gas 
at supersonic velocity. As the radius of this front approaches 
the Stromgren radius 



rs 



\47m^ 



qb 



1/3 



(48) 



(where A'ion is the number of ionizing photons produced by 
the source per second and q:b is the case B recombination 
coefficient), it slows and changes to a D-type front. As it 
does so, a shock front forms and detaches from the ioniza- 
tion front, preceding it into the surrounding gas. The sub- 
sequent expansion phase leads to the formation of a dense 
shell of Hi gas between the shock and the ionization front. 
The thickness of this shell is small compared to the ra dius of 
the H II region. The latter is given approximately by (Yorke 



1986) 



n = rs 



1 + 7^ 

4 rs 



to) 



4/7 



(49) 



Here Cs is the isothermal sound speed in the Hll region, 
and t — to is the time since the formation of the shock. The 
expansion phase ends when the H 11 region attains pressure 
equilibrium with the surrounding gas. This occurs at a ra- 
dius 



n ^ [ — ) rs, 



(50) 



where Ti and T2 axe the temperatures of the Hi and Hll 
gas respectively. For gas cooled by H2, the former will be 
approximately 200K (or smaller if HD cooling is important), 
while the latter will be approximately 10* K; thus rf ~ 20rs. 
However, unless the gas density is extremely large, the main 
sequence lifetime of the central star will generally be too 
short to allow the H 11 region to reach this pressure equilib- 
rium; in this case, the maximum size of the Hll region is 
given by ri, evaluated at the end of the star's life. 

If, rather than a uniform density distribution, we have 
an instead an inhomogeneous distribution, with uniform 
mean density (n) , then the above calculations still give us a 
reasonable estimate for the size of the resulting Hll region, 
as long as we replace by (n^). We generally express this 
by means of a clumping factor C such that C = {n^)/{n)^. 
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This results in a value for rs a factor C^''^ smaller than in 
the uniform density case. 

The other simple case we examine is the propagation of 
an H II region into a radially decreasing density distribution. 
This problem has be en studied by Franco, Tenorio-Tagle & 
Bodenheimer (199C), who show that if the density drops 
off more sharply than n oc r~^^^ then the resulting If ll 
region is unbounded by recombination. This implies that the 
entire halo may become ionized; however, as this requires at 
least one ionizing photon per hydrogen atom (assuming that 
secondary ionization by energetic photons is negligible), this 
cannot take place on a timescale less than 



tic 



M 



3.2 X 10*®M 



(51) 



where A'^ion is the number of ionizing photons produced per 
second and M is the mass of gas in the halo. 

To apply these calculations to our model halos, we need 
to specify A^ion- We can calculate this in a similar fashion to 
A'dis, by modelling the star as a black body. Whilst this will 
inevitably be inaccurate, the values we obtain compare well 
with thos e from the more detailed treatment of Tumlinson 
& ShuU (2000). For the case of a 25 M© star considered 
previously, we find that Mon — 7 x 10** s~^. If we take the 
temperature of the H ll region to be 10'' K, then 



rs — 60n pc, 



(52) 



where we have used the case B recombination coefRcient 



from Osterbrock (1988) 



(53) 



For our uniform density halo model, this implies 

8.1 X 10" 

= , .3 pc. 
(1 + z)^ 

However, given the unrealistic nature of this density profile, 
this figure is probably not relevant. For the more realistic 
truncated isothermal sphere model, we have 



3.8 X lO-" 



pc. 



(54) 



assuming that the Hll region is confined to the core of the 
density profile. From this we can calculate ri, given some 
value for t — to- We choose to adopt a value of 2 Myr, cor- 
responding to a typical OB star lifetimej^ In figure ^ we 
plot ri (in units of the core radius ro) for halos of total mass 
lO'^ M0 and 10* Mq, at a range of redshifts. Also plotted is 
the radius at which the truncated isothermal sphere density 
profile becomes steeper than r'^''^; this occurs at approxi- 
mately 2.9ro. 

For halos in which ri > ro, our analysis, which assumes 
constant density, begins to break down; for ri > 2.9ro, it 
breaks down completely, and the H 11 region will in reality 
become unbounded. We see from figure that this is true 
for all but the most massive halos. Taking this at face value, 
it suggests that, as well as photodissociating H2 throughout 
a halo, a single massive star will also ionize the entire halo. 

t This is an overestimate, but will be reasonably accurate so long 
as the initial R-type phase is much smaller than the subsequent 
D-type phase. 




Figure 10. The final size of the H II region formed by a 25Mq star 
in halos of total mass 10^5 Mq (dotted line) and 10* (dashed 
line) that is predicted if the density is assumed constant and equal 
to the central density of a truncated isothermal sphere of the same 
mass, forming at the specified redshift. It is plotted in units of ro, 
the core radius of the truncated isothermal sphere. Also plotted 
is the radius at which the density profile of the sphere begins to 
fall ofi' more steeply than r~^^^ (solid line). 



However, this prediction may owe much to the overly sim- 
plistic nature of our halo model. In a real halo, several factors 
are likely to combine to reduce our estimated values of ri. 
Significant clumping is likely, and will increase the effective 
density and hence decrease n. Cooling will also increase the 
core density above that predicted by the truncated isother- 
mal sphere model. Finally, non-thermal pressures due to in- 
falling gas or turbulence may also restrict the expansion of 
the Hll region. 

Given these uncertainties, plus the rather ad hoc na- 
ture of the adopted density profile, the final size of the first 
Hll regions is rather unclear. On the one hand, they may 
be comparable in size to the halo; on the other, they may 
prove to more closely resemble the ultracompact H 11 regions 
seen in some local molecular clouds. To decide between these 
scenarios will require a better understanding of the density 
and velocity structure of the first protogalaxies than we cur- 
rently possess, a deficiency that high resolution nume rical 
simulations are beginning to address (Abel et al. 2000). 

Finally, it should be noted that even if the final Hll 
region is large, this will significantly affect our conclusions 
only in the smallest halos. In larger halos, tdis < tion, so 
cooling is still regulated in the first instance by photodisso- 
ciation. Similarly, clumps dense enough to preserve their H2 
until collapse will in general have < tian, and thus will 
only be affected by the Hll region if they lie close to the 
massive star. 



5.4 Dust 

Dust, if present, can affect our simple picture in three ways. 
Firstly, H2 can form by direct association on grain surfaces. 
The rate of H2 formation by this process is independent of 
ionization. Conse quent ly, at low temperatures the argument 
of Nishi & Susa ( 1999 ) is no longer valid, and the H2 abun- 
dance may be significantly larger than 10"''. 

Second, dust absorbs strongly in the Lyman-Werner 
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bands, and will reduce the number of photons able to cause 
photodissociation. 

Finally, dust implies the presence of metals, which may 
replace H2 as the dominant coolant at low temperatures. 
This means that H2 dissociation need not imply cessation of 
star formation, if the abundance of metals and dust is high 
enough to cool the gas. 

With each possible effect, we associate a critical metal 
abundance, above which the effect becomes important. We 
estimate these below. 



5.4-1 H2 formation on grains 



H2 forms on grai n surfaces in the local ISM at a rate (Drains 
&: Bertoldi 19961) 



Rdust = 6.0 X 10 



-18t-,1/2 2 -3 -1 

T ' Uh cm s , 



(55) 



where nn is the H I number density. If we assume both that 
dust in low metallicity gas has the same distribution of grain 
sizes as that in the Milky Way, and that the ratio of dust 
to metals is independent of metallicity, then the rate of H2 
formation on grains in metal-poor gas is simply 



-RduBt = 6.0 X 10"'°r 



18rT-,l/2 



z 



2 -3 - 

tih cm s 



(56) 



Both of these assumptions are somewhat ad hoc but they 
are the simplest that we can make, and should give us a 
reasonable idea of the importance of dust. 

If destruction of H~ by means other than H2 formation 
is negligible, as will be the case prior to star formation, then 
we can directly compare 7?dust with the rate of the H" gas 
phase reaction (equation We see that at low metallicities 
the gas phase reaction will dominate, while at higher metal- 
licities formation on grain surfaces dominates. A reasonable 
definition of the critical metallicity is that at which the two 
rates are equal. This is given by 



2'crit 



Td = O"d,1000-^H ( 



(59) 



where (Td,iooo is the absorption cross-section at 1000 A, Nil 
the column density of hydrogen, again assuming a constant 
dust to metal ratio. Settin g — 1, and adopting a value of 
2 X 10~ cm for crd,iooo 

( Praine fc Bertoldi 199^ ), we can 
solve for the critical metallicity at which the dust becomes 
optically thick. We find that 



Zcrit = 5 X IO'VjT'^G 



(60) 



The critical metallicity thus depends upon the column 
density, and consequently upon the properties of the individ- 
ual halo or clump. The range of column densities included in 
our model is large, being approximately 10^^ to 10^^ cm~^ 
for halos and 10^^ to 10^^ cm for clumps. Adopting the 
median values of 10^^ cm~^ and 10^^ cm~^ as representative 
of 'typical' halos and clumps, this implies critical metallic- 
ities of 0.5 Zq and 5 x 10"'^ Zq respectively. This suggests 
that dust absorption in halos is unlikely to be important un- 
til substantial enrichment has occured, while clumps become 
shielded at much smaller metallicities. 



5.4-3 Alternative coolants 

The presence of metals makes available a large number of 
alternative cooling processes, such as fine structure emission 
from O I and C 11, or rotational emission from CO. At a suffi- 
ciently high metallicity, these dominate the low temperature 
cooling. To properly assess the metallicity at which this oc- 
curs is a complex problem, lying far beyond the scope of 
this paper. However, a study along these lines has recently 
been performed by Omukai (200C). He finds that cooling 
from sources other than H2 typically becomes dominant at 



(57) 



a critical metallicity Zcrit — Zq. While the situation he 
considers - protostellar collapse to very high density, with 
no UV radiation field - differs markedly from our own, the 
differences are unlikely to make alternative coolants signifi- 
cantly more effective. 



For a protogalaxy with temperature T = 10 K and ioniza- 
tion a;o = 3 X 10~* (corresponding to the residual abundance g CONCLUSIONS 
from recombination), this implies that 



1.6 X W'Z, 



0- 



(58) 



Thus, H2 formation on dust may become important at 
metallicities as small as a few thousandths of the solar value. 

Once stars have formed, the H~ pathway becomes far 
less effective as much of it is now photodissociated rather 
than forming H2 , while -Rdust is unaffected. However, by this 
stage the bulk of the H2 has already formed and the 'extra' 
H2 forming on dust for Z ~ Zait will increase tdis only 
slightly. 



5.4.2 Extinction 

To assess the importance of dust extinction, we calculate the 
optical depth due to dust at lOOOA. If we again assume that 
primordial dust is similar to that in the Milky Way, and that 
the dust to metal ratio is constant, then the optical depth 
is given by 



In this paper we have examined the destruction of H2 within 
primordial protogalaxies by ultraviolet radiation from the 
first massive stars to form within them. To do this we have 
constructed a simple analytical model of photodissociation 
which incorporates the essential physics and which should 
give a reasonable order of magnitude estimate of the time 
taken to destroy the bulk of the H2. Using this model, we 
have shown that for diffuse gas in small protogalaxies the 
conclusions arrived at by Omukai & Nishi are broadly cor- 
rect: a single massive star will produce sufficient ultraviolet 
radiation to suppress H2 cooling within its own lifetime. On 
the other hand, in larger protogalaxies the photodissocia- 
tion timescale becomes comparable to the stellar lifetime, 
and the effects on gas cooling are correspondingly smaller. 

We have also examined the effects of photodissociation 
upon dense clumps of gas, such as we expect to find in star 
forming regions. Using a simplified clump model, we find 
that there is a density-dependent critical distance, -Dcrit, be- 
yond which the photodissociation timescale is longer than 
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the clump free-fall timescale. For moderate density clumps, 
Dcrit is comparable to the dimensions of the star-forming 
region, but it decreases sharply as clump density increases. 
Since clumps can continue to collapse only if cooled by H2 , it 
is reasonable to assume that the survival of H2 is a prerequi- 
site for star formation. Thus, stars will form only in clumps 
that lie sufficiently far from any massive star. An alterna- 
tive way of considering this is to suppose that stars will only 
form in clumps above some density threshold, where this is 
chosen such that the corresponding Dcrit is much smaller 
than the dimensions of the protogalaxy. In this case, an ap- 
propriate value for the density threshold is approximately 
10^ cm~^, although this is at best a rough estimate. 

We have also examined a number of factors which may 



compliqaLe this simple pictuie. The must impuitauL piuve 
be ioni zing radiation and the prcocnGC of duot and mctala. — 

The effects of ionizing radiation are sensitive to the 
small-scale density distribution within the halo, and are thus 
hard to assess accurately. However, we have shown that even 
if there is sufficient flux to ionize most or all of the halo, this 
will generally take longer than H2 photodissociation. Con- 
sequently, it should have little effect on our central results. 

Dust and metals can lead to more significant changes 
in the physics, but these begin to take effect only once the 
metallicity of the gas exceeds a few thousandths of the solar 
value. 

Our overall picture of star formation in the first proto- 
galaxies is one in which the formation of dense clumps plays 
a central role. Ultraviolet feedback is effective at suppressing 
cooling in diffuse gas, but has little effect on dense clumps. 
Moreover, while it is likely that supernovae will eventually 
remove all of the gas from the system, we expect this to 
occur after clump collapse. An upper limit to the star for- 
mation efficiency is thus set by the fraction of gas that has 
formed into dense clumps by the time that the first mas- 
sive stars ignite. This may be small, but it is too early to 
conclude this with any real certainty, and thus scope for sig- 
nificant star formation within the first cosmological objects 
still exists. 
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